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Abstract 

Following the recent theory of Lovelock-Brans-Dicke gravity, we continue to investigate the conditions 
to support traversable wormholes by the gravitational effects of spacetime parity and topology, which 
arise from the nonminimal couplings of a background scalar field to the Chem-Pontryagin density and 
the Gauss-Bonnet invariant. The flaring-out condition indicates that a Morris-Thorne-type wormhole 
can be maintained by violating the generalized null energy condition, and thus also breaking down the 
generalized weak, strong, and dominant energy conditions; meanwhile, analyses of the zero-tidal-force 
solution show that the standard null energy condition in general relativity can still be respected by the 
physical matter threading the wormhole. This way, the two topological effects have to dominate over the 
ordinary-matter source of gravity, and the scalar field is preferred to be noncanonical. By treating Brans- 
Dicke gravity as a reduced situation of Lovelock-Brans-Dicke gravity, we also examine the Brans-Dicke 
wormholes and energy conditions. 
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Bonnet invariant; flaring-out condition; generalized and standard energy conditions 


1. Introduction 

A wormhole is a fascinating passage as a shortcut connecting two distant regions in a spacetime or 
bridging two distinct universes. Pioneering investigations of wormholes can date back to the Einstein- 
Rosen bridge [1] in general relativity (GR), and earlier constructions of wormholes, such as those con¬ 
verted from the Kerr-Newman family of black holes, suffer from severe instability against small pertur¬ 
bations and immediate collapse of the throat after formation [2]. 

Modem interest in wormholes are mainly based on the seminal work of Morris and Thome on 
traversable Lorentzian wormholes [3], and the way to convert them into time machines [4]. Morris 
and Thorne firstly designed the metric with the desired structures of a traversable wormhole, and then 
recovered the matter fields through Einstein’s equation. It turns out that the energy-momentum tensor 
has to violate the null energy condition, and thus it needs exotic matter to maintain the wormhole tunnel 
[3]. The standard energy conditions, however, are a cornerstone in many areas in GR, such as the classi¬ 
cal black hole thermodynamics [5, 6]. Thus, much effort has been made to minimize the violation of the 
energy conditions and reduce the encounter of exotic matter at the throat (e.g. [3, 7, 8]). 

The search for promising candidates of exotic matter is not an easy job, and only a small handful 
situations are recognized, such as the quantum Casimir effect and the semiclassical Hawking radiation. 
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while all classical matter fields obey the standard energy conditions. With the development of precision 
cosmology and the discovery of cosmic acceleration, various models of dark energy with exotic equa¬ 
tions of state have been proposed, which provide new possibilities to support wormholes, such as those 
supported by the cosmological constant [9], phantom- or quintom-type energy [10, 11], generalized or 
modified Chaplygin gas [12, 13], and inferacfing dark sectors [14]. 

On fhe ofher hand, as an alfernafive to fhe mysferious dark energy, modified and alfernafive fheories 
of relafivisfic gravify beyond GR have been greafly developed to explain fhe accelerated expansion of 
fhe Universe. The higher order terms or exfra degrees of freedom in fhese fheories yield anfigravify 
effecls, which overtake fhe gravifafional affracfion of ordinary maffer af fhe cosmic scale. Lobo look 
Weyl conformal gravify as an example and suggesfed fhaf modified gravities provide anofher possibilify 
lo supporl fra versatile wormholes [15]: if is fhe generalized energy conditions fhaf are violafed, white fhe 
sfandard energy conditions as in GR may remain valid. To dale, Ihis proposal has been applied lo exacl 
solulions of Morris-Thome-lype wormholes in various modified gravilies, such as fhe melric f(R) [16], 
nonminimal curvafure-malfer coupling [17], Brans-Dicke [19], modified feleparralel [20] , mefric-Palalini 
hybrid/(/?) [21], and Einslein-Gauss-Bonnel gravities [22]. 

In Ihis paper, we will look info Iraversable wormholes and fhe sfandard energy condifions in Lovelock- 
Brans-Dicke gravify [23], which lakes info accounl fhe gravifafional effecls of spacefime partly and lopol- 
ogy by fhe nonminimal couplings of a background scalar field fo fhe Chem-Ponfryagin densily and fhe 
Gauss-Bonnef invarianf. This paper is organized as follows. We firslly review fhe gravify fheory in Sec. 2, 
and derive ils generalized energy condifions in Sec. 3. Then fhe condifions to support Morris-Thorne- 
lype wormholes are investigated in Sec. 4, which are extensively examined by a zero-lidal-force solulion 
in Sec. 5. Also, comparison wilh wormholes in Brans-Dicke gravify is sludied in Sec. 6. Throughoul fhis 
paper, we adopf fhe geomelric convenfions T"^ = T®^^, - S^T"^ • • • and = R^'^j^y wifh 

fhe melric signafure (-, -i- -i- -i-). 


2. Lovelock-Brans-Dicke gravity 

Recenlly we have discussed a new fheory of alternative gravify which has been dubbed as Lovelock- 
Brans-Dicke (LBD) gravify [23]. This fheory is given by fhe action 

*Slbd ^ I + >Sm wilh 

I \ 

0 + bgj - y - 2V(cf>) , 

where (p = 0(x“) is a background scalar field, {a, b) are dimensional coupling conslanls (note: b is hatted 
lo be disfinguished from b - b{r) in Secs. 4, 5 and 6, which is a sfandard denolalion for fhe shape 
funcfion in wormhole physics), ojl denoles fhe dimensionless Lovelock parameter luning fhe kinetics of 
(p{x°), V{(p) refers lo a self-inferacfion pofenlial, and as usual fhe matter aclion is given by fhe matter 
Lagrangian densily via S,n - f d^x In Lq.(l), *RR and g denote fhe Chem-Ponlryagin densily 

and fhe Gauss-Bonnef invarianf, respeclively, 

^ _ pA'*' oabyS 

RR — RapysR - 2^a/3fivR ygR , ^2) 

g -.= R^- AR^yR>^^ + R^avpR^'^"'^ , 

where *Raj 3 y 6 '■= ^^apiivR^^yg is the left dual of Riemann tensor, and represents the totally antisym¬ 
metric Levi-Civita pseudotensor with eoi 23 - = 1 / V^- Note that unlike the other two 


-^LBD - 


Itirr 
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curvature invariants [R, ^], the term *RR in ^lbd is divided by this is because *RR itself already 
serves as a covariant density for 5LBD^ as opposed to ^/^R and therein. 

>5lbd is inspired by the connection between GR and Brans-Dicke gravity, and proposed as the 
Brans-Dicke-type counterpart for the classic Lovelock action in Lovelock’s theorem [24], i.e. .Sl - 
J {r - 2A + -^*RR + + Sf„- *Sl is the most general action made up of algebraic 

curvature invariants that yields second-order field equations in four dimensions, and limits the field equa- 
fion to be Einstein’s equation equipped with a cosmological constant A. The Chern-Pontryagin and the 
Gauss-Bonnet invariants in 5 l do not influence fhe field equation, because *RR and are equal to 

the divergences of their respective topological currents (see Ref. [23] and the relevant references therein); 
instead, the nonminimally 0-coupled covariant densities ^*RR and in the LBD action Eq.(l) 

will have nontrivial contributions to the field equafion. Recall fhaf for fhe fwo invarianfs *RR and Q, 
fhe former is related to the spacetime parity with J *RR proportional to the instanton number of the 
spacetime, while the latter’s integral / dx^^ equates the Euler characteristic of the spacetime. 

Hence, EBD gravity has taken into account the gravitational effects of the spacetime parity and the Euler 
topology. 

The extremized variational derivative SSi^bdISs^'' = 0 yields the gravitational field equafion 

(p + (kpvD - V^Vy) (p + 

/ 1 \ 

- ^ (V^0Vy0 - -g^yVQ,</>V"<^j -I- V{4>)g^y = 

where '■= ^ collecfs fhe confribufions from fhe Chem-Ponfryagin densify wifh nonmini- 

mal coupling fo 0(x“), 

v=^ = 2d^cf> • + 6^y„^V“R^J + Idadficp • (*R“/y + , (4) 

and //py^^ := refers fo fhe effecf of extra degrees of freedom from the nonminimally 

0-coupled Gauss-Bonnet invariant, 

- IR (gpyD - V^Vy) 0 + 4R^ “V„Vy0 + AR, “ 

- 4R,,ucj) - 4g^y • + 4R„^^yV^V>, 

with □ := denoting the covariant d’Alembertian. Compared to the field equations of fhe f{R, Q) 

and f{R,Q,5('m) generalized Gauss-Bonnel gravities wifh generic ^-dependence [25, 26], we have re¬ 
moved fhe algebraic terms in by the Bach-Eanczos identity IRR^y - 4R^°‘Rav - ARa^pyR"^ + 

IRfjapyRy = \Qgiiv- Immediately, the trace of the field equafion (3) is found fo be 

-(pR + ^ + {3 + 2bR) ncp - 4Sr"^Vo.V^0 -f AV{(p) = (6) 

where = 2Rn(p - AR^^^Vp(p, and is always fraceless. 

On fhe ofher hand, for fhe scalar field 0(x"), fhe exfremizafion 6 S\jBT)ld(p = 0 directly leads to the 
kinematical wave equation 

2mLn0 = -Ir+ -^*RR + bg]4>+—'^a<P'^‘^<P + 2V^(p, (7) 

I xFg I ^ 
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with g“^djj(py and := dV{(f))ld(p. Along with the trace equation (6), 

it yields the dynamical wave equation 

[icoL + 3 + IbRjacf) = -+ bgj4> + SnT^’^'> + + 2V^(f> - 4V{4>) , (8) 

which explicitly relates the propagation of 0(x") to the trace of the matter tensor for the energy- 
momentum distribution. 

In this paper, we will work out the conditions to support traversable wormholes in LBD gravity by 
the nontrivial gravitational effects of spacetime parity and topology due to the nonminimal ^-couplings. 
To begin with, we firstly derive the generalized energy conditions for LBD gravity. 


3. Generalized LBD energy conditions 

In a region of a spacetime, for the expansion rate 9(o of a null congruence along its null tangent 
vector field and the expansion rate 6(„) of a timelike congruence along its timelike tangent u^, 6(i) and 
9(u) respectively satisfy the Raychaudhuri equations [6] 

K{e)9(t) - ^ , (9) 

K{uAu) - ^ - o'pvcrj") + - RfivU^u''. (10) 

The inaffinity coefficients are zero K(£-) = 0 = K(u) under affine parameterizations, the twist vanishes 
cOfjyOf'^ = 0 for hypersurface-orthogonal foliations, and being spatial tensors = 0 = cr^^yU^^ the 

shears always satisfy cr^ycr^'' > 0. Thus, to guarantee d9(()ldA < 0 and d9(u)ldT < 0 under all circum¬ 
stances - even in the occasions 9{f) = 0 = 0(„), so that the congruences focus and gravity is always an 
attractive force, the following geometric nonnegativity conditions are expected to hold: 




d9(t) 

dA 

d9(u) 

dr 


RA^r > 0 , > 0 . 


( 11 ) 


Note that although this is the most popular approach to derive Eq.(l 1) for its straightforwardness and 
simplicity, it is not perfect. In general 9(f) and 0(„) are nonzero and one could only obtain ^9^^^+RA^^'' ^ 
0 and -I- R^yU^u'' > 0. Thus, it is only safe to say that Eq.(ll) provides the sufficient rather than 
necessary conditions to ensure d9(f)ldA < 0 and d9(u)ldT < 0. Eortunately, this imperfectness is not 
a disaster and does not negate the conditions in Eq.(ll); for example, one can refer to Ref. [27] for a 
rigorous derivation of the first inequality in Eq.(l 1) from the Virasoro constraint in the worldsheet string 
theory. 

On the other hand, consider generic relativistic gravities with the Eagrangian density Aota\ - 
(R, Rf^yR/^^Ai ,■■■,'&, -I- A,i, where R; ^ % {gafi , Ryj^y/}, '^yRf,av/i ,■■■) refers to a generic cur¬ 

vature invariant beyond the Ricci scalar, and d denotes a scalarial extra degree of freedom unabsorbed 
by An- The field equation reads 


'H^y = with ‘K^y 


1 

v=g dgi^'' 


( 12 ) 


where total-derivative terms should be removed in the derivation of ‘K^y. In the spirit of reconstructing 
an effective dark energy, Eq.(12) can be intrinsically recast into a compact GR form by isolating the Ricci 


4 






(13) 


tensor R^y out of ^^y'. 


Rf^y - Ugf^y ^ ^^G.sTjf^ with - ^G^y - SnGTj^^^ 

^ tJeff 

where Ggff denotes the effective gravitational coupling strength, and it is recognized from the coefficient 
of the matter tensor t‘'™\ refers to the total effective energy-momentum tensor, and = 

- r^y \ with collecting all the modified-gravity nonlinear and higher-order effects. Thus, all 
terms beyond GR have been packed into and Geff. 

Following Eq.(13) along with its trace equation R = -STrGeffT^®^ and the equivalent form R^y = 
S/rGeff the geometric nonnegativity conditions in Eq.(ll) can be translated into the 

generalized null and strong energy conditions (GNEC and GSEC for short) 

G^sTf^^ > 0 (GNEC) , +> 0 (gsec), (14) 

where = 0 for the GNEC, and = -1 in the GSEC for compatibility with the metric signature 
(-, -I- -I- -I-). We further supplement Eq.(14) by the generalized weak energy condition 

GeffT^y^^M^M’' > 0 (GWEC), (15) 

and the generalized dominant energy condition (GDEC) that GeffT^y^M^n'' > 0 with GeffT^y^u^ being a 
causal vector. 

Note that for the common pattern of the field equations in modified gravities, we have chosen to 
adopt Eq.(13) rather than R^y - \Rgiiv - where G is Newton’s constant. That is to say, 

we do not absorb Geff into T^y^ so that GeffT^y*^ = G7^y^; as a consequence, Geff shows up in the 
generalized energy conditions as well. This is because the effective matter-gravity coupling strength Geff 
plays important roles in many physics problems, such as the Wald entropy of black-hole horizons [28] 
and the cosmological gravitational thermodynamics (e.g.[29]), although the meanings and applications 
of Geff have not been fully understood (say the relations between Geff and the weak, Einstein, and strong 
equivalence principles). 

Geff and vary among different theories of modified gravity, which concretize Eqs.(14) and (15) 
into different sets of generalized energy conditions. Eor EBD gravity summarized in Sec. 2, we have 


Geff-<^-^ and ^ 


(16) 


with the components of T^y®^ given by 


%7:T^ify = (V^ Vy - g^yO) ^ ^ - Eg^y . 


Hence, for EBD gravity, the GNEC, GWEC and GSEC are respectively 


0-+ V^Vy0 + y V^,^Vy0 - aH%^> - | > 0, 


/(CP) jirrm) 


(17) 


(18) 
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(p ^u^u'' 




+ v, 






-bH\ 


(GB) 


liv 


+ 4 >- 


□0+^V„0V> +V|>0, (19) 


and |87rr^“^ + V^Vv</. + 

+ ^(p-^(sKT^’^^ + - (1 + 2bR)a((> - 2vj > 0, (20) 

while the GDEC can be concretized in the same way. Among all generalized energy conditions, Eq.(18) 
clearly shows that the GNEC is not influenced by the background potential V = V{(f>) of the scalar field. 

Particularly, EBD gravity reduces to become GR for the situation ^(jc") = G~^ = constant and 
V{<p) = 0, as *RR and in Eovelock’s action .Sl do not affect the field equation. Then Eqs.(14) and 

(15) reduce to become the standard energy conditions for classical matter fields [5]: 

Tj^^er>0 (NEC) , > 0 (WEC) , ^ (SEC). (21) 


4. Conditions to support wormholes in LBD graity 

4.1. Generic conditions supporting static, spherically symmetric wormholes 


If has been nearly three decades since the classical work of Morris and Thome, and nowadays the 
Morris-Thorne metric for static spherically symmetric wormholes is still the most useful and popular 
ansatz to study traversable wormholes. The metric reads [3] 

ds^ = -e^‘^(’-Ut^ + {\- 


b{r) 




( 22 ) 


where (l)(r) and b{r) are the redshift and the shape functions, respectively, and the radial coordinate r > ro 
ranges from a minimum value ro at the wormhole throat to infinity. 0(r) is related to the gravitational 
redshift of an infalling body, and it must be finite everywhere to avoid the behavior q ^nd 

consequently the existence of an event horizon. b{r) determines the shape of the 2-slice {t - constant, 6 - 
7r/2) in the embedding diagram; it satisfies b{r) < r to keep the wormhole Eorentzian, b{rQ) = ro at 
the throat, and b{r)lr ^ 0 at r ^ oo if asymptotically flat. Moreover, the embedding of the 2-slice 
ds^ “ yields the geometrical “flaring-out condition” (b - b'r)lb^ > 0, which 

reduces to become b'{ro) < 1 at the throat r - ro with b{ro) = ro [3]. Here and hereafter the prime 
denotes the derivative with respect to the radial coordinate r. 

Eollowing the metric Eq.(22), in the null tetrad adapted to the spherical symmetry and the null radial 
congruence. 





r J 



m^ - 


V2, 


(0.0.1.^). (23) 


one could find the outgoing expansion rate 0(/) and the ingoing expansion rate 6(„) to be 


dp ) - - (pnp + Pnp) - - 
r 



1 

b { n ) - Pnp + Pnp - - 
r 



(24) 
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where pnp := and pnp := nfm^'^yn^ are two Newman-Penrose spin coefficients. Thus the 

metric ansatz Eq.(22) guarantees that the spacetime is everywhere anti trapped as Gq) - 20(„) > 0, which 
is a characteristic property of traversable wormholes and resembles white holes [30]. Also, Eq.(24) 
shows that the expansion rates are independent of the redshift function d>(r), and the spacetime is free of 
apparent horizons for r > vq. 

Since the outward-flaring constraint (b - b'r)I >0 solely comes from the embedding geometry, it 
is independent of and applicable to all gravity theories. In GR through Einstein’s equation, this condition 
implies that all infalling observers threading a Morris-Thorne wormhole will experience the violation 
of the standard null energy condition > 0 [3, 4]. Similarly, according to the GR form of the 

field equation (13), the flaring-out condition implies that wormholes in EBD gravity are supported by 
the breakdown of the EBD generalized energy conditions as in Eqs.(18) and (19). On the other hand, in 
principle it may still be possible to preserve the standard energy conditions in Eq.(21). Thus, to fulfill 
the constraint {b - b'r)lb^ > 0 in EBD gravity, a possible way to violate the GNEC while keeping the 
standard NEC, i.e. < 0 and > 0, can be 

0 < . (25) 

As another example, violation of the GWEC and preservation of the WEC, i.e. < 0 and 

> 0, can be realized if 

0 < < u^u'' - V^Vv0 - y V^0Vv<^J - + □(/> + eJ . (26) 

Eqs.(25) and (26) indicate that and which represent the effects of the spacetime parity and 

topology, jointly with should dominate over the material source of gravity. Also, a noncanonical 
scalar field (tUL < 0) is preferred than a canonical one (wl > 0) to help support the wormhole. 

Note that in Eqs.(25) and (26) we have assumed 0“' = Ggff > 0. This is inspired by the fact in /(/?) 
gravity that the effective coupling strength Ggff = df(R)ldR =: /r has to satisfy /« > 0 to guarantee 
that in the particle content via the spin projectors, the graviton itself and the induced scalar particle are 
not ghosts [31]. Similarly, in scalar-tensor gravity [f((p)R - h{(l))W(p - 2U{(p)] + ££m in 

the Jordan frame, Geff = f{(pT^ should also be positive definite so that the graviton is not a ghost [32]. 
More generally, for modified gravities of the field equation (13), an assumption Ggff > 0 can not only 
simplify the generalized energy conditions Eqs.(14) and (15), but also help reduce the violation of these 
conditions. 

In Sec. 5, we will demonstrate by a zero-tidal-force solution that Eq.(25) can really be satisfied 
while Eq.(26) is partially falsified for the same numerical setups. To facilitate the discussion, we further 
concretize the tensorial inequalities Eqs.(25) and (26) into an anisotropic perfect fluid form. 

4.2. Supporting conditions in anisotropic fluid scenario 

In accordance with the nonzero and unequal components of the Einstein tensor G^y, one can as¬ 
sume an anisotropic perfect-fluid form T^y - diag [-p(r), P''(r), R^(r), R^(r)j for T^y^^ and each of its 
components. Here T^y is adapted to the metric signature (-, -i- -i- -i-), with p standing for the energy 
density, P’ for the radial pressure, and for the transverse pressure orthogonal to the radial direc¬ 
tion. In wormhole physics, it is P'" that helps to open and maintain the wormhole tunnel, so in the 
context below we will be more concentrative on P’’ rather than P^. Then the generalized energy condi¬ 
tions in Sec. 3 imply Geff(peff + Pg^) > 0 for the GNEC, GeffPeff > 0 and GesipeS + P^g) > 0 for the 


7 


GWEC, GeffCpeff + Pgjf + > 0 and Geff(peff + ^ 0 for the GSEC, as well as GgffPeff ^ 0 and 

GeffPeff > |Geff^eff| for the GDEC, with Ggff removable when Ggff > 0. 

In fact, the perfect-fluid form of T^y clearly shows that violation of the null energy condition - that 
is to say, giving up the dominance of the energy density over the pressure, will imply the simultaneous 
violations of the weak, strong, and dominant energy conditions. This chain of violation happens for 
both the standard and the generalized energy conditions, and in this sense, it is sufficient to consider the 
violation of the null energy condition. According to the GNEC in EBD gravity, it requires + 

^eff^ < 0 to make wormholes flare outward, with pgff = p,„ + + apcp + ^Pgb and + 

aP^p -I- ^PQg; under the Morris-Thome metric, we have pcp = 0 = P[^p, 


%np^ ^ 1 


<!)>' -P f 


+ 


2(l>' mL 

r 2 (p 


+ (b - b'r) + V , 

2r2 


(27) 


SttP 


r 





(28) 


Sn-pGB = -T 


- 4br(p"{b — r) — 2(p'{2r - 3b){b - b'r) + (bb'r + 2r^ — b'r^ + b^ - 3br'j 

+ 20'>V2 {bb'r - b'r^ - b^ + br) + {b - rf + {r^ - Sbr + 4b^) 

+ (t>'(p'r {?ir^ - \6br + 4b'r^ + b'^r^ - 6bb'r + 9b^'j , and (29) 


- ^5 


(p' {b - b'r) {3b - 4r + b'r) + 2^"br{r -b) + 4O'0V (ffiV + 2){b - r)^ 


- Sttpgb ■ (30) 


5. Zero-tidal-force solution 

In this section we will continue to work out an exact solution of Morris-Thorne wormholes in EBD 
gravity, so as to better analyze the flaring-out condition for the wormhole throat, and examine the states 
of the generalized and standard energy conditions. 

There are two functions to be specified in the Morris-Thome metric Eq.(22). To be more concen- 
trative on the wormhole throat and the embedding geometry, we will consider a zero redshift function 
<l'(r) = 0 or e2<i>(P = which corresponds to vanishing tidal force and stationary observers [3]. In this 
situation, the EBD curvature invariants and the Einstein tensor read 

ow 

R - —^ , *RR = 9 = Q, G^y = r~^ ■ diag[ - b'r, -b, b - b'r, b - b'r{, (31) 

and thus the componential field equations G^y = directly illustrate the influences of the 

flaring-out condition {b - b'r)lb^ > 0 to - diag [-peff, PIq, 

To simplify the dynamical wave equation (8), we assume the potential V{^) to satisfy the condition 
-2V, which integrates to yield 


V{(P)^Vocp\ 


(32) 


where Vo is an integration constant. Moreover, we adopt the following power-law ansatz for the static 






and spherically symmetric scalar field, 


<t>{r) = 00 



(33) 


where 0o and the power index A are constants, and vq is the throat radius vq = min(r). 

With these setups, the kinematical wave equation (7) leads to 

(2 + wlA) rb' - a>LA(A - l)b + mLA(A - 2)r - 4Vo0o j = 0. (34) 

Solving this equation for b{r) with the boundary condition b{r - tq) = vq, we obtain the shape function 


2Vo0o?‘n 


b{r) = 


(^) -(^) 


A-3 


oJi^A^ — 2 cl)i^A + a — 3 


+ 


culA (A - 2) r - 2ro (y) 

wlA^ — 26UlA — 2 


^lA( 1 -A) 

iU|^A +2 


and thus 


2Vo0o?‘o 


b - b'r -- 


(> 


_l_ (~) — (A — — 


(A-2)(^) 


A-3 


2ro 


cJi^A{l-A) 

^0\ ‘^L^+2 


colA^ ~ 2ojlA + a — 3 ojlA + 2 \ r 

At the throat r = ro, the flaring-out constraint (b - b'r)lb^ is evaluated as 


4Vo0O?‘o + 

^ - 

Ct>LA + 2 


< 1 


(35) 


(36) 


(37) 


There are five parameters in b{r), among which {ro, 00? A, Tq) attribute to our solution ansatz for the 
homogeneous scalar field 0(r) and the potential T(0), while a»L comes from LBD gravity. To illustrate 
the wormhole geometry, we will adopt the following setups for these parameters. 


(1) Without any loss of generality, let ro = 1 for the throat radius, and 0o = 1. 

(2) According to Eq.(33), asymptotic flatness of the spacetime requires A > 0 so that the scalar field 
monotonically falls off as r ^ oo; moreover, 0(r) is positive definite and meets the expectation 
Geff = 0“^ > 0 for the effective gravitational coupling strength, so that the graviton of LBD gravity 
is non-ghost for the sake of quantum stability. 

(3) A repulsive potential hill L(0) > 0 tends to open and maintain the wormhole tunnel, while a 
trapping potential well L(0) < 0 would collapse the wormhole tunnel. Thus, in our numerical 
modelings, let To - 1 > 0 so that L(0) serves as a potential hill. 


(4) Furthermore, it follows from Eq.(37) that the Lovelock parameter summarily satisfies 

2 

ruL < - — < 0 for ro = 00 = Vo = 1 and A > 0. 


(38) 


This agrees with the indication of Eqs.(25) and (26) that a noncanonical (a>L < 0) scalar field could 
best help support the wormhole. 


With the numerical setups in Eq.(38), only two parameters wl and A remain flexible in determining the 
behaviors of b{r) and {b-b'r)lb^, where A tunes the spatially decaying rate of the scalar field; appropriate 
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values of ojl and A should validate b(r) < r, b - b'r > 0, and wl < -j < 0. In Fig. 1, b{r) is plotted 
at the domain r > ro = I, and the wormhole solution Eq.(35) is confirmed to be Lorentzian. In Fig. 2, 
we plot b - b'r and equivalently verify the flaring-out condition {b - b'r)lb^ > 0. In both figures, we fix 
A = 2.3 and illusfrafe fhe dependence on wl (nofe fhaf inside Figs. 1 ~ 4, a>L is femporarily wriffen as oj 
for fhe sake of greafer clarify). 



r 


Figure 1: With ro = = Vo = 1, A = 2.3 and in the domain r > tq = 1, b(r) is plotted as the solid curves for various oil, along 

with the dotted diagonal for the auxiliary function b(r) = r. For ujl = {-2, -4, -6, -8,-10 • • •) < -2/A = -2/2.3 in light of 
the numerical setups in Eq.(38), b(r) always falls below the auxiliary diagonal b(r) = r. Thus, 1 - b(r)/r is positive definite 
and the wormhole solution Eq.(35) is Lorentzian. Moreover, the curve b(r) approaches the dotted line when £Jl goes to -oo, 
i.e. lim b{r)/r - 1. 


2.5 

2 

1.5 

1 

0.5 

0 

Figure 2: With ro = (/iq = Vq = 1, A = 2.3 and in the domain r > rg = l,b- b'r is plotted for cul = {-2, -4, -6, -8, -10 • • •) < 
-2/A = -2/2.3 and manifests itself to be positive definite. This equivalently confirms the outward-flaring condition (b - 
b'r)lb^ > 0 of the embedding geometry. Moreover, the curve b — b'r tends to coincide with the horizontal r-axis when £Jl 
approaches -oo, i.e. lim b - b'r = 0, which is consistent with the tendency lim b(r)/r = 1 in Fig. 1. 

Wifh fhe Einsfein fensor G^y given by Eq.(31) and (f>~^ > 0, adding up fhe componenfial field equa- 
fions G‘f = -87r0“^peff and G'^^ = one could obfain b - b'r = -r^ ■ 87r0“^(peff + Thus, 

for fhe numerical sefups summarized by Eq.(38), Eig. 2 nof only verifies fhe positive definifeness of 
b -b'r, buf also implies fhe violation of fhe GNEC 0“^(peff + < 0 - and consequenfly fhe GWEC, 

GSEC and GDEC in EBD gravify. On fhe ofher hand, can fhe sfandard energy condifions in Eq.(21) sfill 



r 
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- 0 , 02 -*^-^^^^^^^^^^^- } 

1 1.001 1.002 1.003 1.004 1.005 1.006 

r 


Figure 3: With ro = (^o = Vo = 1, A = 2.3, b = -\ and in the domain r > ro = 1, 8;r(p„, + is plotted as the solid 
curves, while the dotted horizontal depicts the zero reference level. The first subfigure shows that for ujl = {-4, -6, -8,-10), 
8/r(p,„ + P^) is positive definite with the expected asymptote lim 8;r(p,„ + PJ„) = 0^, so the standard NEC is respected. However, 

r—»co 

as oil further decreases, 8;r(p„ + gradually falls below the dotted horizontal near the throat r ^ tq = 1, which has been 
illustrated for oji = {-30, -40, -50, • • • , -1000) in the second subfigure by magnifying the region tq = 1 < r < 1.0006. Thus, 
large negative values of ojl (numerical analysis gives ojl £ -12.9) are unfavored in light of 8;r(p,„ + P'^) > 0. 


hold along the radial direction for the matter threading the wormhole? The energy density p,„ and the 
radial pressure vary for different types of physical matter, and p„j + PJj, relies on the the equation of 
state PJ„ = Thus, we choose to calculate p^ + from an indirect approach. Considering that 

= -ipm +p<f, + apcp + bpGB) and = P;;, + P^ + aP^p + bP[^^, Pm and P^„ can be recovered by 

Skp,„ ^b'rcp- Snp^ - SttBpob , - ^^bPcB ’ (^ 9 ) 

where, according to Eqs.(27)-(30) with <l)(r) = 0, we have 

8w = (i - + V T^) ^ 

^^P^ = 1^" + (41) 

2 

S^PGb - (b - b'r) {3b - 2r) + 2br(f>"{r - b)^ (42) 

and SttP^b = ^{b - b'r){b' - 2). (43) 

In Fig. 3, ^n{pm + P'm) plotted as the solid curves, where we let S = -1 < 0 for the Gauss- 
Bonnet matter-topology coupling strength so that the Gauss-Bonnet part of LBD gravity could yield 
antigravitational effect to help maintain the wormhole tunnel. p,n + P^ is positive definite for wl - 
{-4, -6, -8, -10) < -2/A and thus the standard NEC pm + P^ > 0 is respected by the physical matter, 
despite the violation of the GNEC due to (6“^(peff + Pg^^) < 0; in fact, this has realized the null-energy 
supporting condition of Eq.(25) in an anisotropic perfect fluid form. However, large negafive values of 
wl is unfavored: careful numerical analysis finds fhaf fhe standard NEC becomes slightly violated, i.e. 
Pm + P'm < 0 for ^ -12.9 in the very close vicinity of the wormhole throat. 
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r r 

Figure 4: With ro = = Vo = i,A = 2.3, b = -I and in the domain r > tq = 1, is plotted as the solid curves for Wl = 

{-4, -6, -8,-10). Although p„, > 0 in the distance, one always observes p,„ < 0 near the throat ro = 1, and thus the violation 
of the standard WEC, SEC and DEC. Moreover, as shown in the second subfigure for Wl = {-13, -23, -33, ■ ■ ■ , -1003), the 
intersection point between 8^p„, and the dotted zero reference level moves leftwards when ujl decreases, so the violation of 
Pm ^ 0 gradually reduces. 




Figure 5: With ro = </io = Vo = 1, = -6, A = 2.3 and in the domain r > ro = 1, we plot 8n(p„, + P''^) and 8np„, for different 

Gauss-Bonnet topology-gravity coupling strength, as is given the decreasing series b = {-2,-4,-6,-8,-10). The standard 
NEC always holds with p,„ + > Q. Moreover, the intersection point between 8np„ and the dotted reference level moves 

leftwards when b <Q increases from £ = -10 to -2„ so the violation of p„ > 0 gradually reduces. 
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Figure 6: With ro = i^o = Vo = 1, £Jl = -6, b = -4 and in the domain r > ro = 1, we plot &n(p„ + and 87rp„ for different 
decaying rate of the scalar field, as is given by the increasing series A = {1.3,2.3,3.3,4.3). The standard NEC always holds 
with Pm + P'^ > 0. Moreover, the intersection point between and the dotted reference level moves leftwards when A > 0 
increases from A = 1.3 to 4.3, so the violation of p„, > 0 gradually reduces. 


Validity of the standard weak, strong and dominant energy conditions requires us to check the posi¬ 
tivity of the physical matter density Plotting Snpm for cul - {-4, -6, -8, -10 • • •) in Fig. 4, we find 
p,n < 0 near the wormhole throat, and the violation of p,„ > 0 can be reduced with the decrement of cul 
in the domain wl < -2/A; actually, this has negated the weak-energy supporting condition of Eq.(26) in 
the anisotropic perfect fluid form, despite the validity of the null-energy Eq.(25). As expected, when one 
goes way from the wormhole throat, the normal behaviors p,„ > 0 and lim p^ = 0^ are recovered, and 

r—»oo 

thus the standard WEC becomes valid as p„j -i- PJJ, > 0 for r > rp = 1 in light of Eig. 3. 

Having seen from Eig. 4 that the decrement of the noncanonical a»L could reduce the violation of 
Pm ^ 0, we cannot help but ask are there any other factors that could help protect the standard WEC? 
The answer is yes. In Eigs. 5 and 6, we respectively fix {wl - -6, A = 2.3} and {wl - -6, b = -4} 
to plot [Snipm + P^), STrpmj. It turns out that when the standard NEC is obeyed, i.e. pm + > 0, the 

increment of b (Gauss-Bonnet topology-gravity coupling strength) in the repulsive domain ^ < 0 and the 
increment of A (decaying-rate index of the scalar field) in the domain A > 0 could both help minimize 
the violation of p^ > 0 near the wormhole throat. 

6. Implication: Wormholes in Brans-Dicke gravity 

In the limits a ^ 0 and ^ ^ 0 for the parity-gravity and the topology-gravity coupling coefficients 
in .Slbd> and in the absence of the potential V{(p), EBD gravity reduces to become Brans-Dicke gravity 
with the standard action [33] 

>Sbd = J' ^Pg |(/>P - ^ -I- Sm , (44) 

where m refers to the Brans-Dicke parameter (in distinction with wl for the Eovelock parameter). The 
gravitational field equation d^BD= 0 and the kinematical wave equation d^BD/d0 = 0 are respec¬ 
tively 





+ ^ jV,,0V,0- :rgf,yVa4>^"<p\ ^ 


„(m) 

^ /jy ' 


(45) 
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( 46 ) 


and 2cb - -(pR H—Vq.^V"^ . 

4 > 

With the trace of the field equation -0/?+ ^ VQ.0V“0+3n0 = STrr*^'"^ Eq.(46) leads to the dynamical wave 
equation (2dj + 3)n0 = however, the kinematical equation (46) is preferred so that temporarily 

we need not worry about for the physical matter. 




Figure 7: With = pQ = I and A = 1.3, we plot b(r) and b - b'r as the solid curves. The first subfigure shows that for 
u) = {-1,0,1,2,3 • • •) > -2/A = -2/1.3, b(r) always falls below the dotted diagonal of the auxiliary function b{r) = r, 
and thus guarantees the Lorentzian signature as 1 - b(r}lr > 0. Moreover, the second subfigure verifies b — b'r > 0 for 
tij = {-1,0,1,2,3,4 • • •) > -1/1.3, so the flaring-out condition {b — b'r)lb^ > 0 of the embedding geometry is satisfied. 


For Morris-Thorne wormholes in Brans-Dicke gravity, consider a zero-tidal-force solution <l)(r) = 0, 

( \A 

7 ) (00 > 0, A = constant) of Eq.(33) for the scalar field. Direcfly 
solving Eq.(46) for b{r) with the boundary condition ft(ro) = ro, or just substituting {Vq = 0, tUL oj} 
into Eqs.(35) and (36), we obtain 


b(r) - r + 


d)A{l-A) 

2r-2ro(^)^^ 
mA2 - 2mA - 2 


and 


b - b'r = 


2ro 

ojA+1 



tuA(l-A) 

cbA+2 


(47) 


In light of the flaring-out condition at the wormhole throat, the parameters {A, £ 0 } have to meet the re¬ 
quirement 

^> 0 ) - 1 - < 1 ^ d>A>-2. (48) 

ojA. + 2 

Note that this condition does not conflict with the culA < -2 in Eq.(38): Eq.(48) comes from Eq.(37) 
with {Eo = 0, a»L 1 -^ oj}, while Eq.(38) specifies Eq.(37) by ro = 0o = = E the choices Eq = 0 and 

Eq = 1 (and also Eq = -1, if one would like to check it), i.e. the potential being vanishing, repulsive or 
attractive, lead Eq.(37) to totally different situations. 
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Figure 8: With ro = = 1 and A = 1.3, we plot Sn(p,„ + P'„) for Cj = {-1,0,1,2,3,4) as the solid curves, which always fall 

below the dotted horizontal for the zero reference level. Thus, the standard NEC is always violated. As cj grows from -1 to 4, 
the curve of 87r(p„, + FJjj) gradually moves upward, so in a sense the violation of Sn(p„ + > 0 can be reduced for greater 

values of d) in the domain d) > —IjA. 




Figure 9: With ro = 0o = 1 and A = 1.3, we plot Sttp,,, for d> = {-1,0, 1,2, 3, 4) as the solid curves. These curves stay above 
the zero reference level for d> = jO, 1,2,3,4) and coincide with it for dj - -1; since the curves for tu = jO, 1,2,3,4) are stickily 
close to each other in the first subfigure, we magnify them at 1.2 < r < 1.45 for greater clarity in the second subfigure, which 
shows that from bottom to top or from left to right, the curves correspond to w = 0, ...4 in sequence. Although the energy 
density is nonnegative for a) > -1, the standard WEC, SEC and DEC still fail as 8pm + P'^ <0 for all d) > -21A by Eig. 8. 

Among the three parameters in b{r), {at, A) jointly determine the wormhole structure, while rq acts as 
an auxiliary parameter; for the same reasons indicated in the proceeding section, we inherit the numerical 
setups {rq = 0q = 1,A > 0) to illustrate the Brans-Dicke wormhole Eq.(47), which implies to > -2fA 
from Eq.(48). To start with, the Eorentzian-signature condition r > b{r) and the outward-flaring con- 
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straints b - b'r > 0 are confirmed in Fig. 7. 

Next, let’s check the states of the standard energy conditions in Eq.(21). For the matter threading the 
wormhole, the energy density and radial pressure can be indirectly reconstructed from the field equations 
S/rpm = b'r(f> - Snpip and = -bcp - 8n-P^, where 


8n-p^ = 11 - - 


r 2 (p ] 2r^ 


b'r) , 


(49) 


b 


a 


( \A 

y] and the zero-tidal-force solution Eq.(47), it follows that 


(50) 


87r(pm -t = 


2cbA{cbA^ +A^ + 4A-2)+ 4(A^ +A-\)r 


cbA{A-l) 


2A(mA-tA-tl) r 


(d)A -t 2)(d>A2 - 2mA - 2) 


r^+3 _ 2 ^a - 2) r^+3 ’ 


(51) 


1 r(m -I- 6) - [ 
8np„, = - -I-- 


. o 3&A(A-1)] 
" + -^+ &A+2 


mA(A-l) 

mA+2 


(mA2 - 2mA - 2) 


(52) 


Based on Eqs.(51) and (52), the behaviors of 8 n'(p„, -i- P'^) and 87 rp„, are illustrated in Eigs. 8 and 9 with 
the numerical setups {ro = (po = EA = 1.3 > 0, m > -2/A). Unfortunately, despite the nonnegative 
energy density, the standard null - and thus weak, strong and dominant energy conditions are always 
violated along the radial direction as pm + Pm <0 for r > ro; to make matters slightly better, Eig. 8 
indicates that in a sense such violation could be reduced for greater values of m in the domain m > -2/A. 
Moreover, as shown in Eig. 10 which fixes m = 1 and studies the influences of A instead, we notice that 
even pm > 0 no longer holds throughout r > ro for a spatially quickly decaying (A > 2) scalar field. 

Comparing Eigs. 8 ~ 10 of Brans-Dicke gravity with Eig. 3, one could find that due to the presence 
of the potential hill V{(p) > 0 and the possibly antigravitational Gauss-Bonnet effect, EBD gravity could 
“better” protect the standard NEC when supporting wormholes. The results in this section supplement the 
earlier investigations in Ref. [19]. Moreover, recall that in scalar-tensor theory with the total Eagrangian 
density .if = [f{0)R - h{(p) ■ - 2U{(p)] + if),, in the Jordan frame, we have the following 

conditions for the sake of ghost-freeness and quantum stability [32] : the graviton is non-ghost if /(0) > 0 
(as mentioned before in Sec. 4.1), while the scalar field 0 (a:") itself is non-ghost if 

+mh(4>)>0. (53) 


Eor Brans-Dicke gravity with f{(p) = (p and h{(p) = m/i/i, it requires <;& > 0 and m > -3/2 to be totally 
ghost-free. Thus, the lessons from Eigs. 8 ~ 10 are consistent with the argument of Ref.[32] that in 
scalar-tensor gravity, there exists no static, spherically symmetric wormholes that are both ghost-free 
and obeying the standard NEC. 
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(a) From top to bottom, A =0.1, 0.4, ■ ■ ■, 1.3. 


(b) From top to bottom, A =1.2, 2.2, 3.2, 4.2. 




(c) From top to bottom, A =1.1, 0.9, ■ ■ •, 0.1. 


(d) From top to bottom, A =2, 3, • • ■, 8. 


Figure 10: With ro = 0o = 1 and w = 1, we plot Sn(p„ + P''^) and 8;rp,„ for different values of A which governs the spatially 
decaying rate of the scalar field. Figs. 10a and 10b, with a similar appearance to Fig. 8, show p„, + P'^ <0 for A > 0 so that the 
standard NEC is violated; moreover, with Fig. 10c in an analogous pattern to Fig. 9, one finds p„, > 0 for slow decaying rate 
A = jO.1,0.3, ■ ■ ■ , 1.1). As the most interesting observation. Fig. lOd shows that p„, is no longer positive definite throughout 
r>ro for high decaying rate A > 2. 
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7. Discussion and conclusions 


In Secs. 4 and 5, we have seen that pcp and did not help in supporting Morris-Thome wormholes; 

this is because identically vanishes for all spherically symmetric spacetimes (no matter static or 
dynamical). In fact, it can be directly verified that the spacetime parity will come into effect via nonzero 
in generic axially symmetric spacetimes, say the metric below that generalizes Morris-Thome into 
rotating wormholes [7] : 




b{r, 6) 


dr^ + r^^dO 


i2 ■ 2 

' + Sin 


9 {dt - ojdip)^ ], 


(54) 


where, as in the Kerr or Papapetrou metric, o) = a>(r, 6 ) is the angular velocity dipidt acquired by a test 
particle falling to the point (r, 9) from infinity. 

Also, the wormhole geometry is not only related to the energy-momentum distribution of the physical 
matter through the gravitational field equation, buf also fo fhe propagafion of Ihe scalar field fhrough fhe 
kinemafical wave equafion. Hence, in Secs. 5 and 6, for fhe sake of simplicify, we have chosen fo 
“recover” fhe shape function b{r) and fhus fhe wormhole geomefry from fhe kinemafics of 0(r), i.e. 
Eqs.(7) and (46), while fhe field equations were employed fo analyze fhe energy conditions. When 
seeking for zero-lidal-force solufions wilh a vanishing redshifl function 0(r) = 0, Ibis provides a simpler 
melhod lhan fhaf in Ref. [19] for Brans-Dicke gravity, or Ref. [21] for hybrid mefric-Palafini f{R) gravify 
which is equivalenf fo fhe mixfure of GR and fhe m = -'ill Brans-Dicke gravify; fhey solve for b{r) from 
fhe dynamical wave equafion (i.e. Klein-Gordon equafion) rafher fhan fhe kinematical wave equafion, 
and fhus have fo involve fhe frace of fhe mailer fensor = -p,„ + P'l^+ IPj^ righl from fhe beginning. 
However, when looking for more general solutions wilh 0(r) 0, one should sfill lurn fo fhe melhod in 

Refs. [19] and [21], as if becomes insufiicienf fo defermine fhe Iwo Morris-Thome functions (<l>(r), b(r)} 
from a single kinemafical wave equafion. 

To sum up, in Ihis paper we have investigated fhe condifions fo support Iraversable wormholes in 
LBD gravify. The flaring-oul condifion, which arises from fhe wormholes’ embedding geomefry and 
fhus applies fo all mefric gravities, requires fhe violafion of fhe slandard NEC in GR and fhe GNEC 
in modified gravities. Moreover, fhe breakdown of fhe null energy condifion simulfaneously violales 
fhe weak, slrong and dominanl energy condifions. Wilh Ihese consideralions, we have derived fhe gen¬ 
eralized energy conditions Eqs.(14), (15), (18) and (19) for EBD gravity in fhe form lhal explicilly 
conlains fhe effective gravifalional coupling slrenglh Geff = 0”^- These energy condifions have been 
used fo conslrucf fhe condifions supporling Morris-Thome-lype wormholes, including fhe lensorial ex¬ 
pressions Eqs.(25) and (26), and Iheir anisofropic-perfecl-fluid forms in Sec. 4.2. Moreover, in Sec. 5 
we have oblained an exacl solution of fhe Morris-Thome wormhole wilh a vanishing redshifl funclion 
and fhe shape function Eq.(35), which is supplemented by fhe homogeneous scalar field (/'(x") = 0(r) 
in Eq.(33) and fhe potential E(0) = Vo<p^. Wilh fhe flexible parameters in Eq.(35) for b(r) specified by 
[00 = fQ = Vo = l,A > < -2/A}, we have furlher confirmed fhe Eorenfzian signalure, fhe flaring- 

oul condition, breakdown of fhe GNEC, and validily of fhe slandard NEC. Einally, we also invesfigaled 
zero-lidal-force wormholes in Brans-Dicke gravify, and have shown lhal fhe condition + p;>ois 
nol so well protected as in EBD gravify. 

Nole lhal nalural existence of dark energy becomes effeclive only al scales greafer lhan IMpc [34]. 
Similarly in modified gravlies, fhe higher-order lerms or exlra degrees of freedom are aslrophysically 
recognizable only al galaclic and cosmic levels. Hence, supporting wormholes by dark energy requires 
fo mine and condense dark energy, while supporting wormholes by modified gravity requires unusual 
disfribulions of ordinary mailer. Eor example, in EBD gravify, fhe joinl effecls of 7/^™^ and 

have fo become dominanl over fhe physical mailer source and fhe scalar field is preferred fo be 
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noncanonical. As a closing remark, we have to admit that wormholes in existent studies are mainly 
theoretical exercises and hypothetical objects, and there seems a long way ahead before wormholes can 
be artificially constructed and put to astronomical use. 
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